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j£ ■ Abstract 
Q 

Dense matter in compact stars cools efficiently by neutrino emission via the direct 
Urea processes n — > p + e-\-i? e and p+e — > n+v e . Below the pairing phase transition 
temperature T c these processes are suppressed - at asymptotically low temperatures 
exponentially. We compute the emissivity of the Urea process at one loop for tem- 
peratures T < T c , in the case where the baryons are paired in the x 5o partial wave. 
The Urea process is suppressed linearly, rather than exponentially, within the tem- 
perature range 0.1 < T/T c < 1. The charge-current Cooper pair-breaking process 
contributes up to 50% of the total Urea emissivity in this temperature range. 

o 



"q ■ 1 Introduction and summary 

The direct Urea process (single neutrino /3-decay) was introduced in astro- 
physics by Gamow and Schoenberg in their studies of stellar collapse and 
supernova explosions [1]. In a newly born neutron star, the proton fraction 
drops to a few percent of the net density at the end of the deleptonization 
phase. At such low concentrations of charged particles the energy and mo- 
mentum can not be conserved in a weak charged current decay of a single 
(on-shell) particle because of the large mismatch between the Fermi energies 
of baryons, therefore the direct Urea process is forbidden. The two-baryon pro- 
cesses are not constrained kinematically [2], and the so-called modified Urea 
process n + n — > n + p + e + z/ e is the main source of neutrino emission at low 
proton concentrations. Boguta [3] and Lattimer et al. [4] pointed out that the 
high density matter can sustain large enough proton fraction (> 11% — 13% 
according to ref. [4]) for the Urea process to work. The proton fraction in dense 
matter is controlled by the symmetry energy of the neutron rich nuclear mat- 
ter which is not well constrained. Exotic states of matter open new channels 
of rapid cooling with neutrino emission rates comparable to the direct Urea 
rate [5]. 
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Below T c pair correlations open a gap (A) in the quasiparticle spectrum of 
baryons 1 . At low temperatures, when the width of the energy states accessi- 
ble to quasiparticles due to the temperature smearing of the Fermi-surfaces is 
less than the gap in the quasiparticle spectrum, the neutrino production by 
the direct Urea process is quenched by a factor exp(— A/T), where A is the 
larger of the neutron and proton gaps. At moderate temperatures T < T c , the 
changes in the phase-space occupation caused by the gaps in the quasiparti- 
cle spectra are non-exponential; and the weak decay matrix elements acquire 
coherence pre-factors, since the degrees of freedom (the Bogolyubov quasipar- 
ticles) are a superposition of the particles and holes describing the unpaired 
state. The phase space changes due to the pairing were cast in suppression fac- 
tors multiplying the emissivity of the processes in the normal state in refs. [7]; 
the role of coherence factors has not been studied so far. 

Our aim below is to compute the direct Urea emissivity in the superfluid 
phases of dense nucleonic matter below T c . We start by expressing the neutrino 
emissivity of the direct Urea process in terms of the polarization tensor of mat- 
ter within the real-time Green's functions approach to the quantum neutrino 
transport [8]. (The derivation in Section 2 applies to a larger class of charge 
current reactions involving /3-decay; an example is the modified Urea process). 
The polarization tensor is then computed at one-loop for the isospin asymmet- 
ric nuclear matter where the S'-wave pairing is among the same-isospin quasi- 
particles. Within the range of temperatures 0.1 < T/T c < 1 the suppression of 
the Urea process will turn out to be linear rather than exponential in temper- 
ature. An additional contribution to the Urea process comes from (previously 
ignored) charge-current pair-breaking process. Both the charge-neutral cur- 
rent counterparts of the pair-breaking processes [9,10,11] and the direct Urea 
processes (suppressed by pair-correlations below T c ) were found important in 
numerical simulations of compact star cooling [12,13,14,15,16,17,18,19]. Be- 
low, we shall treat the baryonic matter as non-relativistic Fermi-liquid; one 
should keep in mind, however, that the relativistic mean-field treatments differ 
quantitatively from their non-relativistic counterparts [20,21]. 



2 Propagators and self-energies 

Stellar matter is approximately in equilibrium with respect to the weak pro- 
cesses - the /3-decay and electron capture rates are nearly equal. Thus, we need 

1 Pairing in higher than / = angular momentum states can lead to quasiparticle 
spectra featuring sections of Fermi surfaces where the gap vanishes; here we shall 
treated only I = 0, S'-wave pairing in the 1 Sq partial wave. The S'-wave pairing 
in the 3 Si — 3 D\ channel is suppressed by the mismatch in the Fermi-energies of 
baryons required by the Urea process [6] and can be ignored in the present context. 
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Fig. 1. The diagrammatic Dyson equation for neutrinos which defines the 
self-energies cr > ' < (q). The free and full neutrino propagators are shown by the 
dashed and double dashed lines, the electrons and paired baryons - by the solid 
and double-solid lines; the wavy lines correspond to the gauge P^ + -boson propaga- 
tor. According to Eq. (2), the neutrino Green's function simultaneously propagates 
a neutrino and anti-neutrino hole; the time-reversed Green's function propagates a 
anti-neutrino and a neutrino hole. 

to compute the rate of, say, anti-neutrino production by neutron /3-decay and 
multiply the final result by a factor of 2. We start with the transport equation 
for anti-neutrinos [8] 



d t + d q u v {q)d x ] U{q) = f ^Tr [a<(q)S>(q) - a>(q)S<(q)} , (1) 

where fp(q) is the anti- neutrino Wigner function, Sp ,<: (q) and a > ' < (q) are their 
propagators and self-energies, uj p (q) = \q p \ is the on-mass-shell anti-neutrino 
energy and q = (qo,q) is the four- momentum. The Wigner function, propa- 
gators and self-energies depend on the center-of-mass coordinate X, which is 
implicit in Eq. (1). If the leptons are on-shell and massless, their propagators 
are related to their Wigner functions via the quasiparticle ansatz [8] 



S ^ = J ni{ 6 (* " U M) fx(q) -S(q + u- x (q)) [1 - £(-?)] ), (2) 

where the index A = v, e refers to neutrino {v) and electron (e) and A = P, e 
to anti- neutrino (p) and positron (e). Both neutrinos and electrons obey a 
linear dispersion relation uj\ = since the electron neutrino mass can be 
ignored on energy scales ~ 10 2 keV, while the electrons are ultrarelativistic at 
relevant densities. In addition to the propagator S£ (q) we will need the "time- 
reversed" propagator S^(q) which is obtained from (2) by an interchange of 
the particle and hole distributions f\(q) <-> I — f\(q). The standard model 
neutrino-baryon charged-current interaction is 



V = -^J {H) J {L \ 4 H) = ^(Cy - C A75 )0, Jfa = ^(1 - 75)^, (3) 

where G = GpcosO with Gp being the weak coupling constant and 9 the 
Cabibbo angle (cos# = 0.973), ip and <p are the neutrino and baryon field 
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operators, cy and ca are the weak charged-current vector and axial vector 
coupling constants. The neutrino self-energies 0" > ' < (q r ) can be evaluated per- 
turbatively with respect to the weak interaction. To lowest order in G the 
second-order Born diagram gives (see Fig. 1) 

- ia>><( qi ) =£^(<7i - q 2 + q) iT» L) tS<(q 2 ) ir{£, iU^>(q), (4) 
im 

where J2 q = / d 4 q/(2n) 4 is the phase-space integration symbol, II^ < (g) is the 
baryon polarization tensor and F^l) = 7^(1 — 7s) is the weak interaction 
vertex, which is independent of q at the energy and momentum transfers 
much smaller than the gauge boson mass. Substituting the self-energies and 
the propagators in the collision integral we find for the loss part [the second 
term on the r.-h. side of Eq. (1)] 

S (q 2 - uJ e (q 2 )) [1 - f e (q 2 )\ 5 (q w + uo^q,)) [1 - II< (q), (5) 

where we retained from the collision integral the only relevant term which has 
an electron and anti-neutrino in the final state. The gain term naturally does 
not contribute to the emission rate. 



3 Direct Urea emissivities 

The anti-neutrino emissivity (the energy radiated per unit time) is obtained 
from Eq. (1) 

d f d 3 q r / \ / \ ,„>. 
£s ll# (?W9) ' (6) 

Carrying out the energy integrations in Eq. (5) and using the equilibrium 
identity nS(g) = 2ig£(go)SmIl^(g), where gB{%) is the Bose distribution 
function and n^(g) is the retarded polarization function, we obtain 

7=2) ip4jis4)l A% + 92 - 9) 

5(uj e + uj v -q G )uj 9 (q 2 )g B (q G ) [1 - f e (u e )] A K (<?i, q 2 )^mU^(q), (7) 



v _^2_ r tC_^L_ 
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Fig. 2. The one loop baryon polarization tensor. The baryon and (amputated) 
VF + -boson propagators are shown by the solid and wave lines. The double arrow 
lines refer to the anomalous propagators (ref. [23], pg. 300). 



where K^{q\ 1 q2) = Tr 



The symbol 3m refers to 



- 7 5 ) siWHi - r) 4*\ 

the imaginary part the polarization tensor's resolvent. Since the baryonic com- 
ponent of stellar matter is in thermal equilibrium to a good approximation, it 
is convenient to use the Matsubara Green's functions (ref. [23], pg. 120) 



G a *> (p, r) = -8 a<7 > {T T a pu (T)al a , (0)), 
Faa'(p,r) = (T r a_ pi (r)a pT (0)), 



(T T at 



t a'_ 



■Pi 



(0)), 



a pa {r) and a pa {r) are the creation and de- 



where r is the imaginary time, al 
struction operators, a =f, J, stands for spin and the propagators are diagonal 
in the isospin space. The quasiparticle spectra of Cooper pairs coupled in 
a relative S'-wave state are (the wave-vectors p and k refer to protons and 
neutrons) e p = ^JQ + A~|, e k = + with (omitting the isospin index) 
£ = p 2 /2m + 3?e £(p) — where //, m and S(p) are the chemical poten- 
tials, masses and self-energies of protons and neutrons. Expanding the self- 
energy around the Fermi-momentum pp leads to £ = p 2 /2m* — //*, where 
\i* = — 5te S(pf) + fi, and m/m* = (1 + <99?e £(p)/<9p| p=PF ). The density and 
spin-density response functions are defined in terms of imaginary-time ordered 
products (T T is the ordering symbol) 



n(g,ig )= E / ^e^a-a'(T T a P+gi(T (r)a p , (7 (r)a p ,_ 9/T ,(0)a p , (7 ,(0)),(9) 

where a = 1 for the vector and cr for the axial vector response (cr is the vec- 
tor of Pauli matrices). The one-loop density-density and spin-spin correlation 
functions are given by 



n 



V/A( 



<li i( lo) = \ E [ G (P> ip) G ( k > ik o) T F{p, ip)F\k, ik ) 

P up 



(10) 



where p = (ipo,p), k = (ip + iqo,p + q), the zeroth components of these 
four-vectors are the complex fermionic Matsubara frequencies, and (3 is the 
inverse temperature. The different signs in Eq. (10) reflect the fact that under 
time reversal the a vertex is even for scalar and odd for spinor perturbations. 
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Note that the one-loop approximation above ignores the vertex corrections, 
which can be treated within the Fermi-liquid theory [22]. The summation over 
the Matsubara frequencies in (10) followed by analytical continuation gives 



Y[R ( \ V- / ( U l U k =F UpUkVpVk V 2 V 2 k TUpU k VpV k \ 

\ a; - e p - + id u + e p + e k + id J I 



where u 2 = (1/2) (1 + |£ p |/e p ) and u 2 + v 2 = 1. The first term in Eq. (11) 
describes the quasiparticle scattering, which is the generalization to superfluids 
of the ordinary scattering in the unpaired state; the second one describes 
pair-breaking, which is specific to superfluids. Taking the imaginary part and 
changing to dimensionless variables y = j3u, x = f3^ p and Wi = /3Aj (i = p, n) 
we obtain 



lining, y) = + I$? A {y) 



(12) 



where the (dimensionless) integrals for the vector (V) and axial vector (A) 
couplings are 



4 C/PB (y)=r dx C^ C/PB (x,y,w p )\E^( X ,y,w n )\ 

J — oo 

[f(±y/tf + Z$ - f(^/x 2 + w 2 + y)] 9(1 - |4 
w n ) = ( w ± e p ) ± e p f - Al\}' 1/2 , and 



x 



where E sc ^ PB (x,y, 



(13) 



4 = sgn(u; ± e p )^ y\(co±e p ) 2 -Al\ - Z 2 } 



(14) 



with I 2 = (p 2 + g 2 )/(2m*) — /i*. The axial- vector integrals I^ C ^ PB are obtained 
by replacing in (13) the coherence factors Cy C ^ PB by C^ C ^ PB , where Cyy^ = 
(upUkTVpVk) 2 , and = u 2 k v 2 ±u p u k v p v k . On substituting Eq. (12) in Eq. (7) 
we obtain 



2n j ^ g B {u) [(I™ + I™) + 3^(/i c + I?) 



d 3 qi 



(27r) 3 2w e 



x 



d 3 q 2 



(2tt) 3 2^ 



[1 - / e (cj e )] 5 3 (q 1 + q 2 ~ q)S(u e + u p - u) ujuj e , (15) 
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Fig. 3. The energy distribution of the neutrino radiation d,J(y)/dy for T/T c = 0.9 
(solid lines), 0.5 (dashed lines), 0.3 (dashed dotted lines). For each fixed T/T c the 
lines correspond to (bottom to top) the pair-breaking, scattering and total contri- 
butions. 

where the (unquenched) axial coupling constant g A = 1.26. Since the neutrino 
momentum is much smaller than the electron Fermi-momentum, q ~ q 1 and 
\qi\ ~ p Fe . With this approximation our final result is 



3G 2 m* n m*p Fe 
e,7 = z J J 



J- 



1 r 

6 J-oo 



e J, 



(16) 



PB\ 
V ) 



3g A (I s A c + lD I dzz'Uz-y). 



The second integral can be expressed through the polylogarithmic function 
Li 4 [— exp(y)], though this is not particularly illuminating. Above T c Eq. (16) 
yields the finite-temperature emissivity of the direct Urea process 



_|_ g ^min r 

dy g B (y) In - - e _ {Xinin+y) J dzz 3 f e {z - y), (17) 

where x min = /3(p^ in /2m* - fj,*) and p min = (m*/q)(uj -//£ + //*- q 2 /2m*) 
[here we assume m* ~ m* = m*\. At zero temperature the logarithm in Eq. 
(17) reduces to y8(— x^), the integrals are analytical and one easily recovers 
the zero-temperature result of ref . [4] . 
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Fig. 4. The neutrino emissivity in units of eo versus temperature (solid lines 
A n (0) = A p (0) = 0.5 MeV, dashed-dotted lines A n (0) = 0.5, A p (0) = 2 MeV). 
The scattering, pair-breaking contributions and their sum are shown by dashed and 
dashed-dotted and solid lines. The upper short-dashed line is the extrapolation of 
the rate for unpaired matter to low temperatures, the lower one corresponds to the 
exponential suppression as discussed in the text. 

4 A model calculation 

To set-up an illustrative model, we assume A n (T = 0) = A P (T = 0) = 
0.5 MeV, m*Jm n = m*/m p = 0.7, 9te E(p Fn ) = -10 MeV and 9te £(p Fp ) = 
—160 MeV at the density 0.24 fm -3 , which is above the Urea threshold in 
charge neutral nucleonic matter under /9-equilibrium for the assumed values of 
the baryon self-energies. The neutron self-energy shift is taken more repulsive 
than the one found in non-relativistic calculations with two-body interactions 
to mimic the effect of repulsive three-body forces. The remaining parameters 
of the model are p Fn = 1-8 fm -1 , p Fp = p Fe = 1.07 fm -1 , /x* = 96.3 MeV, 
H* p = 34.2 MeV and /z e = 212.1 MeV. 

Fig. 3 shows the energy distribution of the neutrino radiation dJ(y)/dy for var- 
ious temperatures. The energy distribution is thermal with the maximum at 
uj ~ 3T in agreement with the fact that each anti-neutrino carries on average 
energy equal to 3T in units of Boltzmann constant. For smaller T/T c the max- 
imum shifts to larger energies corresponding to higher "effective temperature" 
and the peak of the distribution is reduced because of pairing; the thermal 
shape of the distribution remains unchanged. While in the limit T/T c < 1 the 
scattering contribution dominates, at lower temperatures the pair-breaking 
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Fig. 5. The log-log plot of the dependence of the neutrino emissivity in units of eo on 
temperature. The scattering, pair-breaking contributions and their sum are marked 
by circles, squares, and diamonds respectively. 

process becomes increasingly important; e.g. at T/T c ~ 0.3 it contributes 
about the half of the scattering contribution. 

Fig. 4 focuses on the temperature dependence of the direct neutrino emissiv- 
ity in the range 0.1 < T/T c < 1. The important feature here is the nearly 
linear dependence of the emissivity on the temperature in the range 0.1 < 
T/T c < 1; the commonly assumed exponential decay - a factor exp(— A/T) 
with A = max(A n , A p ) - underestimates the emissivity. (Similar conclusion 
concerning the suppression of the direct Urea process by pair correlations was 
reached in ref. [7] which treated the scattering contribution to the emissivity). 
The contribution of the pair-breaking processes becomes substantial in the 
low-temperature range 0.1 < T/T c < 0.4 (about the half of the scattering 
contribution at T/T c ~ 0.3). Increasing the value of the proton gap to 2 MeV 
while keeping the neutron gap at its value 0.5 MeV suppresses the emissivity 
of the scattering process, since at a given temperature the phase space acces- 
sible to the excited states is reduced. The pair-breaking processes are almost 
unaffected since they are related to the scattering of particles in and out-of the 
condensate. The temperature dependence of the net emissivity can be crudely 
approximated as e p = 90 e (T/T c — 0.1); this reduces to the standard Urea 
result for T = T c and sufficiently well (for the purpose of cooling simulations) 
reproduces the numerical result in Fig. 3 for T/T c e [0.1; 1]. Alternatively the 
rate of the Urea process in the unpaired matter can be suppressed by a factor 
(10/9)(T/T C — 0.1). Fig 5 illustrates the low-temperature asymptotics of the 
emissivities for A n = A p = 0.5 MeV. For T/T c < 0.1 the contribution of the 
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pair-breaking process to the total emissivity is equal or larger than that of the 
scattering one. Asymptotically, the logarithmic derivative of the emissivity, 
d\oge u /d\ogT, tends to a constant (negative) value which depends only on 
the magnitude(s) of the zero temperature gap(s), as should be the case for 
exponentially suppressed emissivities. 
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